N = 200;
N1 = 300;
D = 100;
L = 50;
W_true = randn(D,10)*randn(10,L);
X_tr = randn(N,D);
Y_tr = X_tr*W_true;
X_ts = randn(N1,D);
Y_ts = X_ts*W_true;

% missRate = 0.8;
% [N,L] = size(Y_tr);
% Y_tr(find((rand(N,L)<missRate))) = -1;

% [out2.top1,out2.top3,out2.obj,out2.time,out2.iter]  = LEML_missing( X_tr, Y_tr, X_ts, Y_ts,50);
[out1.top1,out1.top3,out1.obj,out1.time,out1.k] = R1MP_f( X_tr, Y_tr, X_ts, Y_ts,200);